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Abstract 

Photoproduction events which have two or more jets have been studied in the W^ p 
range 135 GeV < W 1P < 280 GeV with the ZEUS detector at HERA. A class of events 
is observed with little hadronic activity between the jets. The jets are separated by 
pseudorapidity intervals (At/) of up to four units and have transverse energies greater 
than 6 GeV. A gap is defined as the absence between the jets of particles with transverse 
energy greater than 300 MeV. The fraction of events containing a gap is measured as a 
function of Ar/. It decreases exponentially as expected for processes in which colour is 
exchanged between the jets, up to a value of A 77 ~ 3, then reaches a constant value of 
about 0.1. The excess above the exponential fall-off can be interpreted as evidence for 
hard diffractive scattering via a strongly interacting colour singlet object. 
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1 Introduction 



In high energy hadronic collisions, the dominant mechanism for jet production is described by 
a hard scatter between partons in the incoming hadrons via a quark or gluon propagator. This 
propagator carries colour charge. Since colour confinement requires that the final state contain 
only colour singlet objects, the exchange of colour quantum numbers in the hard process means 
that a jet at some later stage generally exchanges colour with another jet or beam remnant 
widely separated from it in rapidity. Such jets are said to be "colour connected" and this 
leads to the production of particles throughout the rapidity region between the jets. However, 
if the hard scattering were mediated by the exchange of a colour singlet propagator in the 
t-channel, each jet would be colour connected only to the beam remnant closest in rapidity 
and the rapidity region between the jets would contain few final-state particles JTJ. The colour 
singlet propagator could be an electroweak gauge boson or a strongly interacting object, and 
the soft gluon emission pattern produced in each case is similar p|. However the rates could be 
very different. In order to determine the rate of colour singlet exchange processes it has been 
proposed H to study the multiplicity distribution in pseudorapidityQ (r/) and azimuth (<p) of 
the final state particles in dijet events, and to count events with an absence of particles (i.e. 
with a rapidity gap) between the two jets. 

DO @ and CDF || have reported the results of searches in pp collisions at y's = 1.8 TeV for 
dijet events containing a rapidity gap between the two highest transverse energy (E^) jets. 
Both collaborations see an excess of gap events over the expectations from colour exchange 
processes. DO report an excess of 0.0107 ±0.0010 (stat )to^oll(sys. .), whereas CDF measure the 
fraction to be 0.0086 ± 0.0012. We report here the results of a similar search in 7p interactions 
obtained from e + p collisions at HERA. 

In leading order, two processes are responsible for jet production in jp interactions at HERA. In 
the first case, the direct contribution, the photon interacts directly with a parton in the proton. 
In the second case, the resolved contribution, the photon first fluctuates into a hadronic state 
which acts as a source of partons which then scatter off partons in the proton. Fig. 1(a) 
shows schematically an example of colour singlet exchange in resolved photoproduction in 
which a parton in the photon scatters from a parton in the proton, via i-channel exchange of a 
colour singlet object. An example of the more common colour non-singlet exchange mechanism 
is shown in Fig. 1(b). For high E 3 ^ dijet production, the magnitude of the square of the 
four-momentum (|t|) transferred by the colour singlet object is large. Thus it is possible to 
calculate in perturbative QCD the cross section for the exchange of a strongly interacting 
colour singlet object P, 0, [7|, For instance, the ratio of the two-gluon colour singlet exchange 
cross section to the single gluon exchange cross section has been estimated to be about 0.1 ||. 
Studies of rapidity gaps at high \t\ ("hard diffractive scattering") are complementary to studies 
of diffractive hard scattering where the rapidity gap is between a colourless beam remnant, 
produced with low four-momentum transfer with respect to one of the beam particles, and 
hadronic activity in the central detector ||. 

The event morphology for the process of Fig. 1(a) is illustrated in Fig. 1(c). There are two 
jets in the final state, shown as circles in (77, ip) space. Here Ar/ is defined as the distance in 77 

1 rj = — ln(tan ^) where 1) is the polar angle with respect to the z axis, which in the ZEUS coordinate system 
is defined to be the proton direction. 
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between the centres of the two jet cones. For the colour singlet exchange process of Fig. 1(a), 
radiation into the region (labelled "gap") between the jet cones is suppressed, giving rise to the 
rapidity gap signature. Multiplicity fluctuations in colour non-singlet exchange events can also 
produce gaps between jets. In order to disentangle the different mechanisms for gap production 
it is useful to introduce the concept of the 'gap-fraction'. 

The gap-fraction, /(At/), is defined as the ratio of the number of dijet events at this At/ which 
have a rapidity gap between the jets to the total number of dijet events at this At/. For 
colour non-singlet exchange, the gap-fraction is expectated to fall exponentially with increasing 
At/. This exponential behaviour can be taken as a definition of non-diffractive processes Q. 
The expectation follows from the assumption that the probability density for radiation of a 
particle is constant across the rapidity interval between the jets and it is consistent with the 
results of analytic QCD calculations @, and with Monte Carlo simulation (see subsequent 
sections). For colour singlet exchange, the gap-fraction is not expected to depend strongly 
upon At/ H [7|. Therefore, at sufficiently large Arj, such a colour singlet contribution will 
dominate the behaviour of the gap- fraction. The situation is illustrated in Fig. 1 (d), where 
the colour non-singlet contribution is shown as an exponential fall-off, and the colour singlet 
contribution is shown as independent of At]. 

In this paper the gap fraction is studied for a sample of photoproduction events with two jets 
of E^ 1 > 6 GeV. The events are obtained from an integrated luminosity of 2.6 pb _1 of e + p 
collisions measured by the ZEUS detector and have jp centre-of-mass energies in the range 
135 GeV < W w < 280 GeV. Dijet cross sections are measured as a function of Arj for events 
with a gap and for events with no gap requirement. 



2 Experimental Setup 

Details of the ZEUS detector have been described elsewhere [[T(J . The primary components used 
in this analysis are the central calorimeter and the central tracking detectors. The uranium- 



scintillator calorimeter [11] covers about 99.7% of the total solid angle and is subdivided into 
electromagnetic and hadronic sections with respective typical cell sizes, of 5 x 20 cm 2 (10 x 20 
cm 2 in the rear calorimeter, i.e. the positron direction) and 20 x 20 cm 2 . The central tracking 



system consists of a vertex detector [|T^] and a central tracking chamber |13[ enclosed in a 1.43 T 
solenoidal magnetic field. 

A photon lead-scintillator calorimeter is used to measure the luminosity via the positron-proton 
Bremsstrahlung process. This calorimeter is installed inside the HERA tunnel and subtends a 



small angle in the positron beam direction from the interaction vertex |L4] . Low angle scattered 
positrons are detected in a similar lead-scintillator calorimeter. 

In 1994 HERA provided 820 GeV protons and 27.5 GeV positrons colliding in 153 bunches. 
Additional unpaired positron and proton bunches circulated to allow monitoring of background 
from beam-gas interactions. 
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3 Data Selection 



The ZEUS data acquisition uses a three level trigger system. At the first level events are selected 
which were triggered on a coincidence of a regional or transverse energy sum in the calorimeter 
with a track coming from the interaction region measured in the central tracking chamber. At 
the second level a cut was made on the total transverse energy, and cuts on calorimeter energies 
and timing were used to suppress events caused by interactions between the proton beam and 



residual gas in the beam pipe ||15|| . At the third level, tracking cuts were made to reject events 
arising from proton beam-gas interactions and cosmic ray events. Also at the third level, jets 
were found from the calorimeter cell energies and positions using a fast cone algorithm and 
events were required to have at least two jets. 

Charged current events are rejected by a cut on the missing transverse momentum measured 
in the calorimeter. Events with a scattered positron candidate in the calorimeter are rejected. 
This restricts the range of the photon virtuality to P 2 < 4 GeV 2 , and results in a median P 2 of 
~ 10~ 3 GeV 2 . A cut of 0.15 < y < 0.7 is applied on the fraction of the positron's momentum 



which is carried by the photon, where y is reconstructed using the Jacquet-Blondel method |16 
This cut restricts the 7p centre-of-mass energies to lie in the range 135 GeV < W 1P < 280 GeV. 

To select the final jet sample, a cone algorithm [0 is applied to the calorimeter cells. Cells 



within a radius R = y5ri 2 ell + 5ip 2 cell of 1.0 from the jet centre are included in the jet where 
Srj cel1 amd S(j) cel1 represent respectively the difference in pseudorapidity and azimuthal angle (in 
radians) between the centre of the cell and the jet axis. Events are then required to have at 
least two jets found in the uranium calorimeter with E^ 1 > 5 GeV and if et < 2.5. In addition 
the two highest transverse energy jet sQ were required to have Ar] > 2 (i.e. cones not overlapping 
in if) and boost |(r/i + ?7 2 ) I / 2 = |r/| < 0.75. These conditions constrain the jets to lie within the 
kinematic region where the detector and event simulations are best understood. 

To identify gap events, the particle multiplicity is determined by grouping calorimeter cells into 
"islands". This is done by assigning to every cell a pointer to its highest energy neighbour. A 
cell which has no highest energy neighbour is a local maximum. An island is formed for each 
local maximum which includes all of the cells that point to it. The events with no islands of 
transverse energy E^ land > 250 MeV and t] between the edges of the jet cones (as defined by 
the cone radius R) are called gap events. 

A total of 8393 dijet events were selected, of which 3186 are gap events. The non-e + p collision 
background was estimated using the number of events associated with unpaired bunch crossings. 
The beam gas background was found to be less than 0.1%. The cosmic ray contamination is 
estimated to be about 0.1%. For those events in which the low angle scattered positron is 
detected in lead-scintillator calorimeter, P 2 < 0.02 GeV 2 . The fraction of these events is 
around 20%, in agreement with the Monte Carlo expectation. The 43 gap events which have 
At] > 3.5 were also scanned visually to search for contamination from events where the energy 
deposits of the scattered positron or a prompt photon might mimic a jet. No such events were 
found. 



2 In the jets are ordered in pseudorapidity rather than transverse energy and the two jets at lowest 
and highest pseudorapidity are used in the calculation. When the uncorrected gap-fraction is made with this 
selection, it is about 0.01 lower. 
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4 Results 



In section 4.1 we present results obtained from ZEUS data which are not corrected for detector 
effects, by comparing the data to Monte Carlo generated events which have been passed through 
a detailed simulation of the ZEUS detector and selection criteria. The PYTHIA flTSU Monte 



Carlo program was used with the minimum px of the hard scatter set to 2.5 GeV. The GRV [IS 



parton distributions were used for the photon and the MRSA PDA parton distributions were used 



for the proton. Two Monte Carlo event samples were generated. For the first sample ( "PYTHIA 
non-singlet"), resolved and direct photon interactions were generated separately and combined 
according to the cross sections determined by PYTHIA. No electroweak exchange (quark quark 
scattering via ■y/Z or W exchange) events were included. For the second sample ("PYTHIA 
mixed"), 10% of electroweak exchange events were included. This fraction is two orders of 
magnitude higher than the level actually expected from the cross section for these events and 
is chosen in order to mimic the effect of strongly interacting colour singlet exchange processes 
which are not included in PYTHIA. 

In section 4.2 we present the ZEUS data after corrections for all detector acceptance and 
resolution effects. These hadron-level measurements are then compared to model predictions, 
and to the expectation of an exponential suppression of gap production for non-diffractive 
processes. 



4.1 Uncorrected Results 

The energy flow l/NdE^ 11 /d5rj cel1 with respect to the jet axis for cells within one radian in (p 
of the jet axis is shown in Fig. 2(a) for the two highest transverse energy jets of each event. 
PYTHIA mixed events are shown as the solid line. Here and throughout Fig. 2 the data are 
shown as black dots and the errors shown are statistical only. This jet profile shows highly 
collimated jets in the data and a pedestal of less than 1 GeV of transverse energy per unit 
pseudorapidity outside the jet cone radius of 1.0. The pedestal transverse energy is higher 
toward the proton direction. The superposition of profiles of one jet at high rf et and one at 
low rf et leads to the bump at Sr] cel1 ~ 1.5, due to the forward edge of the calorimeter. The 
profiles for the PYTHIA non-singlet sample are not shown, as they are similar to those of the 
mixed sample. The PYTHIA events generally describe the data well, although they are slightly 
more collimated and underestimate the forward jet pedestal. This small discrepancy may be 
related to secondary interactions between the photon remnant and the proton remnant, which 
are not simulated in these PYTHIA samples. Including any kind of multiple interactions in the 
simulation increases the energy flow and particle multiplicity |21| and thus can only decrease 
the number of gaps predicted by the Monte Carlo program. 

The distribution of the total number of events (without any demand on the presence or absence 
of a gap) as a function of At] is shown in Fig. 2(b). It decreases with increasing Ar], extending 
out to Arj ~ 4. The PYTHIA distributions are normalized to the number of events in the data. 
Both PYTHIA samples provide an adequate description of this distribution although the total 
number of events seen at large Ar] is slightly underestimated. 



4 



The distribution of the gap events as a function of Ar/ is plotted in Fig. 2(c) where the 
normalisation for the PYTHIA distribution is the same as in Fig. 2(b). The number of events 
in the data exhibiting a gap falls steeply with Ar/. However the expectation from the PYTHIA 
non-singlet sample falls more steeply than the data, significantly underestimating the number 
of gap events at large Ar/. The PYTHIA sample with a mixture of 10% of electroweak boson 
exchange can account for the number of gap events in the data at large Ar/. However this sample 
significantly overestimates the number of gap events at low Ar/. As mentioned previously, 
including secondary interactions in the simulation could reduce the predicted number of gap 
events and possibly account for this discrepancy. 

By taking the ratio of Fig. 2(c) to Fig. 2(b), the gap-fraction shown in Fig. 2(d) is obtained. 
The gap-fraction falls exponentially out to Ar/ ~ 3.2. Thereafter it levels off at a value of 
roughly 0.08. The PYTHIA non-singlet sample fails to describe the flat region in the data, 
falling approximately exponentially over the whole measured range of Ar/. This sample also 
overestimates the fraction of gap events at low Ar/. The PYTHIA mixed sample can describe 
the flat region of the data but again overestimates the gap-fraction at low Ar/. The gap-fraction 
for the electroweak exchange events alone exceeds 0.4 over the full Ar/ range (not shown). 

The uncorrected data exhibit a flat region at large Ar/ consistent with a colour singlet contri- 
bution of around 10%. Detector effects are expected to largely cancel in the gap-fraction. In 
the next section we find that this is indeed the case and provide quantitative estimates of both 
the discrepancy between PYTHIA and the data and of the significance of the deviation of the 
measured gap-fraction from an exponential fall. 



4.2 Corrected Results 

In order to investigate whether the observed flat region in the gap-fraction might be a detector 
effect, the PYTHIA mixed sample has been used to correct the data for all detector effects, 
including acceptance, smearing and the shift in the measurement of energies. Cross sections 
are determined and the gap-fraction is measured in four bins of Ar/ in the range 2 < Ar/ < 4. 

The cross section da/dArj is measured for dijet photoproduction, ep — > cyp — > eX, where X 
contains at least two jets of final state particles. The cross section is measured in the range 
0.2 < y < 0.85 for photon virtualities P 2 < 4 GeV 2 . The two jets are defined by a cone 
algorithm with a cone radius of 1.0 in (77, </?) and satisfy E J T et > 6 GeV, rj jet < 2.5. The two 
jets of highest EjT f satisfy Ar/ > 2 and \fj\ < 0.75. The rear rf et distribution falls to zero at 
rf et ~ —2, well within the rear detector acceptance. Therefore no explicit rear pseudorapidity 
cut is made. The gap cross section, da gap /dArj, is measured, in the same kinematic range, for 
events with no final state particles with transverse energy E^ arttde > 300 MeV between the 
jet cones. The corrected gap-fraction /(Ar/) is then obtained from the ratio of da gap /dArj to 
da/dArj. 

The efficiency of the data selection described in section 3 for finding events in this kinematic 
region was estimated using the Monte Carlo samples. The combined efficiency of the online 
triggers is at least 80% in every bin of Ar/. The efficiency of the offline selection is about 50% 
leading to a combined efficiency of the online and offline selection criteria of about 40%. The 
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low efficiency of the offline selection is due to the finite detector resolution of the jet energy and 
angular variables, and the steeply falling Ej? spectrum. The shifts and resolutions of these 
variables are consistent with those obtained in extensive studies of the 1993 dijet sample fl22 



The Ej? rUde resolution is 27% with a shift of -14%. The rf art%cle resolution is 0.01 with negligible 
shift. " 

The final correction factors for the inclusive cross section are smoothly varying between 1.6 in 
the lowest A 77 bin and 1.4 in the highest At] bin. The correction factors for the gap cross 
section are between 1.5 and 1.8. The ratios of these correction factors form effective correction 
factors for the gap-fraction which are all within 27% of unity. 

The systematic uncertainties have been estimated by varying the cuts made on the recon- 
structed quantities. The island algorithm for counting particles was replaced by an algorithm 
which clusters cells based on proximity in (rj, (p) space. Also the results of the bin-by-bin cor- 
rection were checked using an unfolding procedure |p3 |. The correction was also performed by 



using the PYTHIA non-singlet sample and by leaving out the leading order direct contribution. 
The uncertainty due to the parton distribution was included. The uncertainty due to a 3.3% 
systematic error in the luminosity measurement was included. Finally the systematic uncer- 
tainty arising from a 5% uncertainty in the mean energies measured by the calorimeter was 
estimated. This represents the largest uncertainty in the two cross sections but cancels in the 
gap-fraction. The largest systematic uncertainty which remains in the gap-fraction comes from 
the variation of the E^ land cut from 200 to 300 MeV. The combined effect of these uncertainties 
is included in the outer error bars in Fig. 3. 

The inclusive and gap cross sections and the corrected gap-fraction as a function of At] are 
presented in Figs. 3(a) to (c) respectively (black dots) and compared with the expectations from 
the PYTHIA non-singlet exchange cross sections (open circles). For the data, the inner error 
bars show the statistical errors and the outer error bars display the systematic uncertainties, 
added in quadrature. The cross section points are plotted at the centres of the bins. The 
gap-fraction points are plotted at the mean Arj values of the inclusive cross section. Numerical 
values for the inclusive cross section, the gap cross section and the corrected gap-fraction are 
provided in Tables 1, 2 and 3, respectively. 

The inclusive cross section is around 5 nb per unit Arj at Arj = 2, falling to about 0.5 nb for 
At] > 3.5. The gap cross section is around 3 nb per unit Arj at Arj = 2 and falls to about 
0.06 nb for At] between 3.5 and 4. The overall normalization of PYTHIA agrees with the data 
within the errors. PYTHIA also describes the shape of the inclusive cross section. However it 
fails to describe the gap cross section, falling too steeply with At] and disagreeing significantly 
in the last bin. 

The corrected gap-fraction falls exponentially in the first three bins but the height of the 
fourth bin is consistent with the height of the third. The height of the fourth bin is 0.11 ± 
0. 02 (stat.)t_Q^l(sys.), which is also consistent with the flat region at large Arj seen in the 
uncorrected gap fraction and inconsistent with the expectation from PYTHIA. 
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5 Discussion 



Two methods have been used to estimate the significance of the excess of the gap-fraction over 
the expectation from multiplicity fluctuations in non-singlet exchange. 

The first method is to take the difference between the data and the PYTHIA non-singlet gap- 
fractions, shown in Fig. 3(c). An excess of 0.07 ± 0.03 is obtained, based entirely on the last 
bin. However this is a model-dependent estimate. For instance replacing the Lund symmetric 
fragmentation function by the Field-Feynman fragmentation function yields a lower predicted 
gap-fraction and a larger excess. Introducing multiple interactions into PYTHIA also lowers the 
fraction of gap events expected. On the other hand, lowering a PT (which controls the hadron 
momentum distribution transverse to the parent parton) in the Monte Carlo simulation from 
the default value of 0.36 GeV to 0.25 GeV produces a gap-fraction that is very like the data. It 
has a height in the fourth bin of 0.07 ±0.02 and therefore if one believes this model, there is no 
significant excess. However this option yields jet-profiles which are narrower than the default 
PYTHIA profiles, which are already slightly narrower than the data. 

The second way to estimate the excess of the gap-fraction over that expected from purely non- 
singlet exchange does not rely on comparisons to Monte Carlo predictions. In Fig. 3(d) the 
gap-fraction is shown again and compared with the result of a two parameter (a,/3) X 2_ fit to 
the expression ff lt (a,/3; Arf) = C(a, j3)e aAv + [3 where C(a,/3) is the normalization coefficient 
constraining ff lt (a,/3; Arj) to 1.0 at Ai] = 2. The result of this fit is shown as the solid curve 
in Fig 3(d), and the exponential (dotted line) and constant (dashed line) terms are also shown. 
The quality of this fit, as indicated by the \ 2 value of 1.2 for the two degrees of freedom is 
superior to that of a fit to an exponential alone which yields \ 2 = 9. The fit parameters 
are a = -2.7 ± 0.3(stat.) ± 0.1(sys.) and (3 = 0.07 ± 0.02(stat.)±° ° l(sys.). The parameter (3 
gives an estimate of the gap-fraction for colour singlet processes. This method uses the full 
information of the four measured data points and is not dependent on the details of the Monte 
Carlo fragmentation model. However, the assumption that the colour singlet gap-fraction is 
constant with Ar) is only one of many possibilities. 

Both the comparison with the default PYTHIA non-singlet prediction and the fit to an expo- 
nential form give an excess of about 0.07 in the gap-fraction over the expectation from colour 
non-singlet exchange. 

The excess in the gap-fraction over the expectation from non-singlet exchange may be inter- 
preted as evidence for the exchange of a colour singlet object. In fact the fraction of events due 
to colour singlet exchange, f(Ar]), may be even higher than the measured excess. As previously 
mentioned, secondary interactions of the photon and proton remnant jets could fill in the gap. 
A survival probability, V, has been defined || which represents the probability that a secondary 
interaction does not occur. Then f(Ar]) = f(Ai]) ■ V. Estimates of the survival probability for 



pp collisions at the Tevatron range from about 5% to 30% |3], |24], ^5|. The survival probability 
at HERA could be higher due to the lower centre-of-mass energy, the fact that one remnant jet 
comes from a photon rather than a proton and the fact that the mean fraction of the photon 
energy participating in the jet production in these events is highP]. Therefore the ZEUS result 

3 The average fraction of the photon energy participating in the production of the two jets |2^| is 0.7 for these 
events. Nevertheless, according to the PYTHIA simulation the dominant contribution in this kinematic regime 
is from leading order resolved events. 
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of ~ 0.07 and the DO and CDF results of 0.0107 ±0.0010( y stat.)t 2 1 l(sys.) and 0.0086 ±0.0012 
could arise from the same underlying process. 

The magnitude of the squared four-momentum transfer across the rapidity gap as calculated 
from the jets is large (\t\ > (£jf*) 2 ). Thus the colour singlet exchange is unambiguously "hard". 

The PYTHIA generator predicts that the ratio of the electroweak (a EW ) to QCD (a® CD ) 
exchange cross sections in this kinematic range is a EW /a® CD < 7 ■ 10~ 4 (compatible with the 
estimation (a/a s ) 2 ). Therefore quark quark scattering via •y/Z and exchange cannot 
explain the height of the flat region in the gap-fraction. On the other hand, using the simple 
two-gluon model for pomeron exchange gives f(Arj) ~ 0.1 ||. Thus pomeron exchange could 
account for the data. 

In summary, dijet photoproduction events with > 6 GeV contain an excess of events with 
a rapidity gap between the two jets over the expectations of colour exchange processes. This 
excess is observed as a flat region in the gap-fraction at large rapidity separation (Ar]= 3.7) 
at a level of 0.11 ± 0.02(stat.)^^(sys.). It can be interpreted as evidence of hard diffractive 
scattering via a strongly interacting colour singlet object. 
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Ar] 


da / dArj 


Statistical 


Systematic 




(nb) 


Uncertainty (nb) 


Uncertainty (nb) 


2.25 


4.93 


0.24 


+0.83 
-0.68 


2.75 


3.06 


0.15 


+0.54 
-0.52 


3.25 


1.67 


0.07 


+0.31 
-0.19 


3.75 


0.54 


0.03 


+0.08 
-0.03 



Table 1: da/dAr] for ep — > cyp — * eX in the kinematic range 0.2 < y < 0.8 and P 2 < 4 GeV 2 
and where X contains two or more jets of E^ et > 6 GeV, rf et < 2.5, \fj\ < 0.75 and Ar] > 2. 



Ar] 


da^P/dAr] 


Statistical 


Systematic 




(nb) 


Uncertainty (nb) 


Uncertainty (nb) 


2.25 


2.85 


0.17 


+0.45 
-0.45 


2.75 


0.66 


0.06 


+0.11 
-0.15 


3.25 


0.16 


0.02 


+0.03 
-0.04 


3.75 


0.06 


0.01 


+0.01 
-0.01 



Table 2: da gap /dAr] for ep — > e^p — * eX in the kinematic range 0.2 < y < 0.8 and P 2 < 4 GeV 2 
and where X contains two or more jets of E 3 T et > 6 GeV, r] jet < 2.5, \fj\ < 0.75 and Ar] > 2 
with no final state particles of Ej? rUde > 300 MeV between the jets. 



Ar] 


/(At?) 


Statistical 
Uncertainty 


Systematic 
Uncertainty 


2.23 


0.58 


0.04 


+0.04 
-0.02 


2.73 


0.22 


0.02 


+0.02 
-0.02 


3.22 


0.10 


0.01 


+0.01 
-0.02 


3.70 


0.11 


0.02 


+0.01 
-0.02 



Table 3: The gap-fraction, /(Ar/), for ep — > e^yp — * eX in the kinematic range 0.2 < y < 0.8 
and P 2 < 4 GeV 2 and where X contains two or more jets of E^ et > 6 GeV, rf et < 2.5, \fj\ < 0.75 
and Ar] > 2. 
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(c) (d) 

Figure 1: Resolved photoproduction via (a) colour singlet exchange and (b) colour non-singlet 
exchange. The rapidity gap event morphology is shown in (c) where black dots represent final 
state hadrons and the boundary illustrates the limit of the ZEUS acceptance. Two jets of 
radius R are shown, which are back to back in azimuth and separated by a pseudorapidity 
interval Arj. An expectation for the behaviour of the gap fraction is shown in (d) (solid line). 
The non-singlet contribution is shown as the dotted line and the colour singlet contribution as 
the dashed line. 
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2 2.5 3 3.5 4 2 2.5 3 3.5 4 

(c) At? (d) At? 

Figure 2: Uncorrected data compared with the predictions from PYTHIA events which have 
been passed through a detailed simulation of the ZEUS detector and of the sample selection 
criteria. The errors shown are statistical only. The transverse energy flow with respect to the 
jet axis is shown in (a) where the data are shown as black dots and the PYTHIA non-singlet 
sample is shown as a solid line. In (b), (c) and (d) the data are again shown as black dots. The 
PYTHIA non-singlet sample is shown as open circles and the PYTHIA mixed sample (which 
contains 10% of colour singlet exchange events) is shown as stars. The number of events versus 
At] is shown in (b). The number of gap events versus At] is shown in (c) and the gap-fraction 
is shown in (d). In (d) the points are drawn at the mean Arj of the inclusive distribution in the 
corresponding bin. 
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Figure 3: ZEUS data (black circles) corrected for detector effects. The inner error bars represent 
the statistical errors from the data and Monte Carlo samples, and the outer error bars include 
the systematic uncertainty, added in quadrature. In (a), (b) and (c) the PYTHIA prediction 
for non-singlet exchange events is shown as open circles. The inclusive cross section is shown in 
(a). The cross section for gap events is shown in (b) and the gap-fraction is shown in (c). The 
gap-fraction is redisplayed in (d) and compared with the result of a fit to an exponential plus 
a constant. In (c) and (d) the points are drawn at the mean A77 of the inclusive distribution in 
the corresponding bin. 
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